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Based on the entropy concept, we define a new clustering order parameter c feasible for finite 
systems of interacting oscillators. Unlike the generalized synchronization order parameters of the 
Kuramoto type, this new order parameter singles out the splay state as the unique state with c = 0, 
thus yielding a positive value whenever there is some kind of cluster formation in the system. It is 
therefore proposed to be monitored alongside the Kuramoto order parameters as a means to quantify 
the overall amount of clustering in the system. 


I. INTRODUCTION 


The phenomenon of collective synchronization among 
coupled oscillators has, and continues to be, an area of 
intense study. Assuming a finite system consisting of 
N oscillators with phase variables 9j, Kuramoto turned 
the intuitive idea of synchronization into the following 
mathematically tractable formula: 
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entropy based clustering order parameter, denoted by c, 
that could be monitored alongside the synchronization 
order parameters of the Kuramoto type. In the pursuit 
of such an order parameter c we will require the following 
properties: 

1. It should mimic r in the sense that it yields 
c = 0 for the splay state and take the value c = 1 only 
for states where all oscillators have the same position on 
the unit circle. 

In addition to this 


The order parameter r thus defined has earned a special 
status since it forms an integral part of the celebrated 
Kuramoto model [ij: 
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2. It should single out the splay state as the unique 
state with c = 0. 

3. It should be devoid of arbitrary parameters (such as 
coarse-graining) and easy to implement numerically for 
a finite system of oscillators. 


which, using the order parameters defined above, can ex¬ 
pressed in the form 

9i{t) = iOt +Krsm[^|; - 9i{t)]. (3) 


However, if the coupling function were to be replaced by 
a higher order harmonic, for example sm[2{9j{t) — 9i{t))], 
then r would no longer adequately quantify the amount of 
clustering in the system. This has called for a generalized 
set of order parameters defined as follows: 
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see Refs. As an example, given a state with two 

equally sized clusters separated by an angle tt we obtain 
r = 0 which is the same value as that obtained for 
the splay state where all the oscillators are uniformly 
distributed over the circle, see Fig. [T] On the other 
hand, for the two-cluster state we obtain r 2 = 1 but we 
still have r 2 = 0 for the splay state. It is clear, however, 
that no Tk will simultaneously single out the splay state 
as the unique state with the lowest order parameter and 
the complete in-phase state as the unique state with the 
highest order parameter. Therefore, in order to achieve 
this objective, in this brief paper we will propose an 


In the infinite-N Kuramoto model d, dl, where the 
oscillators are represented by a density function p there 
is a functional that clearly distinguishes the uniformly 
incoherent state po{9) = ^ from all other states, namely 
the entropy 


p\n[p)d9, (5) 

see Refs. [3,[1|- Arguably, in the finite-N model the state 
corresponding to po ought to be the splay state. However, 
if we were to try to adapt these concepts directly to the 
finite-N model we would face the problem of how to define 
the density of oscillators at any particular angle. This 
would necessitate coarse-graining of the unit circle using 
some mesh size 69 with the density defined accordingly 
as p = 5N/59, with 5N being the number of oscillators 
residing in the interval 69. Obviously, in order to produce 
meaningful outcomes this mesh size must not be chosen 
too big nor too small, moreover, no matter how small we 
make it will never be able to single out the splay state 
nor the complete in-phase state as unique. Faced with 
these problems, we have found a solution by applying the 
entropy concept instead to the set of angle separations 
between the oscillators as will be described in SectionlTTl 



2 


oooo 


r — 1 
r2 = 1 
c = 1 


r = 0 
r2 = 1 
c = 0.61 


r = 0 
r2 = 0 
c = 0.39 


r = 0 
r2 = 0 
c = 0 


FIG. 1. The value of the order parameters r, r 2 and c (see Eqn. O for various cluster formations in a system consisting of 6 
oscillators. The leftmost figure shows the complete in-phase state for which we obtain r = r 2 = c = 1. The rightmost figure 
shows the splay state for which we obtain r = r 2 = c = 0. Among these order parameters, c is the only one that distinguishes 
all the different cluster formations and shows a gradual decline as we move from the complete in-phase state to the splay state 
via the intermediate two-cluster and three-cluster states. 


II. DEFINITION AND ELEMENTARY 
PROPERTIES 

Starting at an arbitrary point on the unit circle we can 
order the oscillators as we encounter them while moving 
along the circle in the positive direction, see Fig. [2] If 
two or more oscillators happen to be at precisely the 
same angle they could either be given an arbitrary order 
among themselves or simply be considered as one single 
oscillator (though still remembering the fixed system size 
N). We let A9 denote the angular separation between 
two consecutive oscillators ordered this way. Then we 
have that 


V A6» = 27r. 
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Furthermore, we define the entropy S of the set of angle 
separations as follows: 
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where 
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In the case A9 = 0 we let P{A9) take its limiting value 
^0) = 0. By applying the Gibbs inequality (see Ref. 
[8|) one can show that S attains its maximum when A9 
is constant, which for a system of N oscillators means 
that A9 = 2tt/N. The latter is the definition of the 
splay state, for which we obtain S' = 1. Moreover, since 
= /3(27r) = 0 and in the interval 0 < A9 < 2-k we 
have that P{A9) > 0, is is clear that the state with all 
oscillators having the same position on the unit circle is 
the only state for which we obtain S = 0. Thus, we have 
arrived at a definition of c satisfying all the requirements 
imposed on it in Section [71 

(9) 



FIG. 2. Illustration of the set of angle separations {A0}. We 
let A6 denote the angular separation between two consecutive 
oscillators as we encounter them while moving along the circle 
in the positive direction. 


In Fig. [T]we have calculated the value of c, alongside r 
and r 2 , for various cluster formations in a system of 6 
oscillators. 


III. SIMULATIONS ON THE KURAMOTO 
MODEL 

In Fig. Owe present the time averaged values of r and 
c as a function of K resulting from computer simulations 
performed on the original Kuramoto model 


Ut) 


-- UJi 


K 


N 


^sm[9j{t) - 9,{t)]. 
i=i 


with N = 100 and the natural frequencies uji drawn from 
a Gaussian distribution with standard deviation 0.1. The 
length of each simulation was set to 10000. As we can see, 
in this case c behaves qualitatively similar to r. However, 
if we modify the coupling function to consider instead the 
system of equations 
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9,{t) = 0 ;,-f — ^sin[2(6»j(t)-6»,(t))], (11) 


1=1 


c= 1-S' 
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FIG. 3. Time averaged values of the order parameters c and 
r as a function of K obtained from simulations performed on 
the original Kuramoto model, Eqns. (1101) . In this case, the 
qualitative behavior of c is similar to that of r. 




FIG. 4. Time averaged values of the order parameters c and 
r as a function of K obtained from simulations performed on 
the modified Kuramoto model given by Eqns. dl). In this 
case c still shows an increasing level of clustering as a function 
of K whereas r does not show any such tendency. 


then, as evidenced by Fig. |4j in this modified system 
c still shows an increasing level of clustering as a func¬ 
tion of K whereas r does not show any such tendency. 
Continuing the comparison one might ask what time av¬ 
eraged value of c we expect for sub-critical coupling, for 
example when K = 0. In the case of r one concludes that 
r ^ 1/^/N for KT = 0 which can be understood through 


the following heuristic argument. The formula for can 

be expressed as: 
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If we assume no coupling or correlation between the 
phases and simply draw them randomly from the uni¬ 
form distribution 9j G t7[0, 27r] then we immediately ob¬ 
tain the expectation value 

(r^) = l/N. (14) 

Attempting a similar analysis on c would bring us far out¬ 
side the scope of the present paper. Nevertheless, by per¬ 
forming a Monte-Carlo simulation we see that for uncor¬ 
related 9j G 17[0, 27r] and sufficiently large N > 100 the 
dependence of c on iV follows closely the curve aj In(A^), 
with a ~ 0.424, see Fig [SI These values can also be 
confirmed by solving the actual differential equations on 
systems of uncoupled oscillators. 



FIG. 5. The red asterisks show the mean value of c as a 
function of N obtained from Monte-Garlo simulations using 
uncorrelated Qj G t/[0, 27r]. The blue curve shows the function 
Q:/ln(A), with a = 0.424. 
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